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We study the propagation of sound-like perturbations created by a jet moving with supersonic 
velocity through the quark-gluon-plasma created in heavy-ion reactions within the model MACE 
(MAch Cone Evohition). Predictions for heavy-ion reactions at RHIC energies (Au+Au-collisions) 
and for Pb+Pb reactions at the LHC (y/s = 5.5A TeV) are presented and potential observations by 
the STAR, PHENIX and ALICE experiments are discussed. 
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I. INTRODUCTION 

First experimental data at the Relativistic Heavy Ion 
Collider (RHIC) have shown the disappearance of back- 
ward correlations in hard processes at central Au+Au- 
Collisions for high-n^-particles [3, Q. Correlation data 
with lower p^-threshold for associated partióles show a 
broader peak or even a double-peaked structure in the 
backward región 0, H i, i, 0, H 0. Three-particle- 



correlations 1C| seem to support the theory of conical 
emission [3, 11, [Hj]. Recently, the observations have 
been confirmed by preliminary data from the CERES- 
experiment at the Super Protón Synchrotron (SPS) [l3l |. 

Conical structures in heavy ion collisions have been 
postulated first by Scheid and Greiner already in 
1974 [l4| and later also by Hofmann, Baumgardt, Stócker 
and collaborators in the second half of the 1970s [HI, [TU, 
13, El, E M EL SI SI SI Hl. In that time, strong 
shock waves in cold nuclear matter were discussed. New 
interest has been rising in this field in recent years, trig- 
gered by the above-named experimental results and new 
theoretical work from Stócker and Casalderrey-Solana et 
in 2004 Ü6|, [13 , and lots of further considerations in 
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^,|4J,|43|,|44j,l45|,|46l. 
The common idea is that an ultra-relativistic jet which 
is moving through an extended Quark-Gluon-Plasma will 
cause some kind of perturbation, maybe sound-like, and 
thereby lose its energy. Then, it cannot be seen in the 
experiments anymore. Since the speed of sound is consid- 
crably smaller than the speed of light, the perturbations, 
if indeed sound-like, should interfere constructively along 
a cone, the so-called Mach cone. 



The term "backward hemisphere" is meant relative to 
one high-energetic trigger-jet. The jet that creates the 
conical structure is supposed to be opposite to the trig- 
ger, i.e. in exact backward direction. 

The opening angle ¡5 of this cone is in static médium 
given by 



sin 



es 



(1) 



where es is the speed of sound in the médium and v is 
the velocity of the jet. The cone is supposed to directly 
emit (or indirectly enhance emission) of partióles. Those 
partióles will be found in a circle around the direction of 
the jet at 
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In two-particle correlations this scenario will mani- 
fest itself as a double-peaked structure in the backward- 
hemisphere. 

The case of non-static médium has been studied in 
several publications as well for homogeneous, non-static 
médium as for simplified of full 3+1-dimensional hydro- 
dynamical evolutions of the underlying médium. All of 
the latter works have described the sound waves and 
Mach cones within the framework of hydrodynamics and, 
for that reason, together with the underlying médium. 

In the presented work, we develop a model that prop- 
agates sound waves through the velocity field as created 
by applying the hydrodynamical equations and discovers 
Mach cone structures automatically. Backreaction of jet 
and sound waves onto the médium is neglected. 



II. HYDRODYNAMICS 
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A model widely used to describe heavy-ion reactions 
is fluid- or hydrodynamics. It has been predicted as a 
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key mechanism for the creation of hot and dense matter 
very early [l5|, [l8[ . Using this approach one assumes that 
the matter described is, at least nearly, in local thermal 
cquilibrium. This is only justified after a short period 
of time. After that, a locally cquilibrated system may 
be formed. The part of the reaction before equilibration 
cannot be described by usual (1-ñuid-) hydrodynamics. 

Therefore, an additional model has to be applied for 
the creation of the first equilibrated state, the so-called 
initial state. This can in principie be any kind of non- 
equilibrium model. The creation of the first cquilibrated 
state can takc different times, depending on the mech- 
anism with which it is reached. Times of the order of 
to ~ 1 fm are reasonable. 

When the initial state is defined, hydrodynamics start 
to work. Depending on number and kind of the assumed 
symmetries in the initial state the fluid development may 
be solvable analytically or only numerically. Within this 
stage of the calculations, one has to assume an equation 
of state (EoS). It must also be chosen if the evolution de- 
scribes a perfect fluid that is perfectly equilibrated at any 
point or if small deviations are endorsed. In the latter 
case, several additional parameters like the heat conduc- 
tivity, shear- and bulk-viscosity have to be introduced. 

In the case of a perfect liquid, there are six equations 
which have to be solved. Among them are one equation 
for the conservation of the baryon number density and 
four for the conservation of energy and momentum: 



d»N»= (3) 

8^= , (4) 

where the interesting variables are given by 

= gu» (5) 

= (e + P)u^u v - Pg» v . (6) 



with the local rest frame baryon density g, the local rest 
frame energy density e, the pressure P and the flow ve- 
locity u^. 

The sixth equation needed is the EoS which gives a 
connection between pressure and the densities: 

P = P(e,g) . (7) 

Within the hydrodynamical framework the propaga- 
tion of sound waves can be described. The speed of sound 
cg can here be derived from the EoS as 

4-£ ■ 

III. INGREDIENTS 

The model developed for this work is independent of 
the model used for the hydrodynamical evolution. This is 
because the hydro-evolution is only used as an input for 
the model, not an inherent part of it. Indeed, any model 



that produces a velocity field can be used. However, of 
coursc, a different input will alter the results. 

Here, we use the Partióle in Cell Code (PiC) which 
bases on the method with the same ñame develo ped 
by Amsden and Harlow in the early 1960s ¡47|, HH 
and upgraded to ultra-relativistic energies by Nix and 
Strottman in the 80s and 90s ¡49l IHol Í5l{ . The method 
takes advantage of a combination of Eulerian and La- 
grangian solution methods; the pressure is calculated 
on an Eulerian grid (where the cells are fixed in coor- 
dínate space), while all current transfers are calculated 
in a larger number of Lagrangian cells (where the cells 
comoving with the fluid) . 

The initial state for the hydro-evolution is calculated 
in an effective ID Yang-Mills coherent field model, where 
the transverse plañe is split up to "streaks" like in the 
"Firestreak" -model of the 1970s [52l[53j. During the evo- 
lution, an ideal gas EoS is used, so that the pressure P is 
connected to the energy density e in a very simple way: 

P = 1/3e , (9) 

and the speed of sound is given as 

cg = 1/V3ps 0.577 . (10) 

This corresponds to a first-order approximation onto 
a Quark-Gluon-Plasma (QGP). The fact that there is 
no phase transition is justified because the main part of 
sound wave evolution will take place in this first, hot and 
dense stage of the reaction. 

IV. MACE 

The model "MAch Cones Evolution" (MACE) is built 
to model the propagation of sound waves through a re- 
alistic, relativistic médium. Our approach is to take an 
existing hydro-evolution and impose a jet and the waves 
it creates as perturbations. Therefore, no backreaction 
of jet or sound waves onto the médium can be modelled. 

A. Jet 

Only the away-side-jet (the jet in the backward hemi- 
sphere) is considered, since the near-side-jet is not af- 
fected by and does not affect the médium. 

The jet starts at the first timestep calculated by the 
hydrodynamical code. This is, depending on the system 
considered, after few fm after the first collisions. 

The jet is created at a given point with a probability 
that is proportional to the energy-density at that point. 
The jet's direction is totally random. The only cut that 
is applied is the exclusión of all jets that will not end up 
in the detector. That is, for RHIC-data we take all jets 
with \r¡\ < 0.7 (corresponding to the STAR-TPC), and 
for LHC-data we cut at \rj\ = 0.9, which is the acceptance 
of the ALICE-TPC. 
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The jet is considered to be high energetic enough not 
to change its velocity and direction. Calculation of jet 
quenching is not within the scope of the model. It henee 
propagates in a straight line with speed of light through 
the médium. At every timestep, it excites sound waves at 
its position. These waves are not enhanced in direction of 
the jet, but are undirected. The direction of the jet will 
be used for correlation considerations in the end. When 
out of the médium, the jet does not excite sound waves 
anymore. 



B. Sound waves 

No premature assumptions on the shape of a result- 
ing mach front can be made in an unforeseeable, inho- 
mogeneous and non-statical médium. We will therefore 
propágate the single elementary waves and identify the 
wave fronts independent of the propagation. This allows 
for parts of the elementary waves to become a part of a 
wave front only after some time and possible deflection. 
Also, some parts may be swamped away out of the wave 
front. 

So, in order to model the sound waves we créate a lot 
of logical partióles, so-called wave markers, at the posi- 
tion of the jet. The word "logical" refers to the fact that 
we do not assign any physical quantities to these mark- 
ers yet; they only represent the position and direction 
of the wave. Random directions will be assigned to the 
markers, so that after a short propagation in homoge- 
neous médium one elementary wave should indeed be a 
spherical wave. Such an elementary wave is created at 
each timestep at the position of the jet. Between the 
timesteps all "waves", i.e. all wave markers, are propa- 
gated. 

The propagation of the wave markers is straightfor- 
ward: The only assumption made is that they move with 
the speed of sound relative to the fluid wherever they are. 
Their direction is adjusted by relativistically adding their 
initial velocity, v, to the flow-velocity of the underlying 
médium at the current point, ü. This initial velocity 
v is calculated by multiplying the direction of the wave 
marker v with the speed of sound: v — csv. 

When a marker crosses the border of a fluid cell, then 
its current propagation fi — r + v ■ dt will be finished 
with the oíd velocity, i.e. no deflection happens at the 
border of two cells. This is justified because the timestep 
dt is much smaller than the dimensions of the cells da:, 
dy and dz. 

If a marker leaves the system, it is deleted. No particle 
emission at this point is assumed, and the wave is also not 
reflected. Fig. [T] shows the position of markers after 10 
timesteps of an exemplified propagation (here, the data 
for a peripheral collision at RHIC have been applicd). 

It should be noted that, by conserving the number of 
wave markers per elementary wave, their "amplitude" 
will decrease with the 1 / r 2 -behaviour expected for spher- 
ical waves. 
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FIG. 1: (Color Online) Wave markers for a jet propagating 
through a peripheral RHIC-collision after 10 timesteps. The 
jet was started near the center of the collision with a rapidity 
of y = 1.1. The red line shows the trajectory of the jet, the 
blue crosses and green stars show the wave markers started 
at odd and even timesteps, respectively. 

Since the underlying hydro-evolution is calculated us- 
ing an EoS with a constant speed of sound (see Eq. ÍTÜty ). 
it would be inconsistent to use a density-dependent speed 
of sound for the propagation of the sound waves. There- 
fore, propagation always happens with es = l/v3. 



C. Mach cones 

In order to see the collective phenomenon of a Mach 
cone, the región oceupied by wave markers is reshaped 
with a set of lines, so-called wave lines. They go along the 
surface of that región, starting at the position of the jet. 
A handy analogy to this are longitudes going from the 
pole of the earth along its surface and thereby reshaping 
it. 

The wave lines are equally distributed azimuthally 
around the jet's direction. Their nodes are positions of 
marker partióles. By restricting the possible nodes to 
those marker partióles that are in the correct azimuthal 
slice the azimuthal position of the wave lines is fixed. 

The sheer surface of the área might be very distorted 
or irregular due to numerical reasons. For instance, due 
to the finite density of timesteps, the surface will contain 
large portions of the same elementary waves (in the limit 
of infinitely many elementary waves each of these only 
contribute one point on a circle in a plañe orthogonal to 
the jet's direction). Also, the random direction of the 
wave markers and their finite number may cause distor- 
tion. It might be the case that in the área where a given 
elementary wave should represent the surface there is no 
wave marker to determine a possible node. 
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FIG. 2: (Color Online) Sketch on how the signáis are created 
from the wave lines. The red dotted line shows the jet trajec- 
tory, the blue dash-dotted line a wave line. The arrows are 
perpendicular to the wave line in constant intervals along it. 

In MACE, the surface is smoothed to elimínate these 
effects. Care is taken not to conceal physical effects dur- 
ing this procedure. 

When all nodes are collected and the lines are com- 
plete, particle flow is considered to be orthogonal to the 
wave lines. The lines are now taken as is. To find the 
places on which to insert a signal we go along the wave 
lines, starting from the jet, and insert a signal after fixed 
distances. Momentum is inserted in the plañe where the 
wave line is, orthogonal to the jet direction (see Fig. [2]). 

The signáis are taken directly and mapped into corre- 
lation functions. There is no new velocity or momentum 
profile created using signal and underlying hydrodynam- 
ical fields. Therefore, background-subtracted data are 
obtained more easlily. The correlation plots shown in 
section Ivl are normalizad to the number of signáis. 

Obviously, the position and density of wave markers is 
totally lost during this procedure. Instead, the density of 
wave lines is important. It is easy to see that this drops 
with the inverse of the distance to the jet axis. Since the 
amplitude of a Mach cone decreases with the same law, 
no additional, artificial decreasing of the signal has to be 
performed. Therefore, it is justified to keep the distance 
between to signáis along the wave lines constant. 

The región right behind the jet is treated as any other 
región in the model. Henee, the cones really have a sharp 
top. Non-linear behaviour and a rounder shape in this 
región is not included. 

V. RESULTS 

The correlation functions obtained with different veloc- 
ity fields created from different initial states show a very 
common shape. Both in central RHIC- (see Fig. [3]) and 
central LHC-Collisions (see Fig. 0} (-y/s^ = 130 GeV 
and y/sÑÑ = 5.5 TeV, respectively) the averaged cor- 
relation functions show a clear double-peaked structure 
in the backward hemisphere. Also, in all cases those 
two peaks are not at the expected positions Atp = ti ± 



FIG. 3: (Color Online) Two-particle correlations (back- 
ward hemisphere) for central gold-on-gold RHIC-collisions 
(V S NN ~ GeV and b = 0). The plot shows the aver- 
age over 1 000 arbitrary jets. The arrows show the "perfect 
mach positions" Aip — iv ± cos^ 1 cg. 
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FIG. 4: (Color Online) Two- (a) and Three- (b) particle corre- 
lations (backward hemisphere) for central lead-on-lead LHC- 
collisions (^/sjsfisf = 5.5 TeV and 6 = 0). Both plots show 
the average over 1 000 arbitrary jets. The arrows show the 
"perfect mach positions" Aip — ty ± eos" 1 cg. 



eos -1 cg, which would be Aípi = 2.2 and Aip 2 =4.1. In- 
stead, they are systematically at positions further apart, 
so the angle to the backward direction is systematically 
bigger. The magnitude of this shift is also very constant 
among the considered cases; it is about SAip w 0.2. Cor- 
relations for mid-central RHIC collisions (not shown) do 
not show a significantly different picture. 

Three-particle correlations (see Fig. [4] (b)) show the 
expected four-peaked structure, but also here it can be 
seen that these peaks are not at the expected position. 

A trigger on the direction of the jet (see Fig. [5]) gives 
the same picture. Both the midrapidity jets and the for- 
ward jets result in correlation functions with peaks at the 
same positions as the peaks from the untriggered data. 
Please note that since all plots are made for central colli- 
sions, in which the system is azimuthally symmetric, an 
additional trigger on in-/out-of-plane-jets is meaningless. 

While triggering on the rapidity of a jet is experi- 
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FIG. 5: (Color Online) Two-particle correlations (back- 
ward hemisphere) for central LHC-collisions for forward- and 
midrapidity-jets. Both lines show the average over 250 jets 
from arbitrary positions. The arrows show the "perfect mach 
positions" Aip = 7r ± eos -1 cg. 




FIG. 6: (Color Online) (a): Two-particle correlations 
(backward hemisphere) for central LHC-collisions for single 
midrapidity-jets started from different positions in the colli- 
sion zone. (b): The starting positions for the jets shown in 
(a). 



mentally a very easy task, it is impossible to model- 
independently measure the origin of a jet inside the reac- 
tion zone. But since all data with averaged origins show 
the same behaviour, it is interesting to look for cases 
where the origin is fixed. 

In Fig.[ü](a) the correlation functions for three different 
jet origins (which are shown in Fig. [ü] (b)) are shown. 
All jets go in midrapidity and into the same azimuthal 
direction. The latter should not make any difference due 
to the symmetry present in central collisions. 



The figure shows how much the correlation functions 
differ. While the correlation function from the centrally 
started jet is both symmetric around ir and has the peaks 
at the expected positions Aip w n ± 0.96, the more pe- 
riphcrally started jets have one of the peaks much smallcr 
than the other and far away from the expected position. 
When averaging over lots of jets coming from periphcral 
positions, it is obvious why the correlation function looks 
like the ones shown in Figs [3l [4] and O 

One important result is that no correlation function 
shows the "wide-peak"-structure as observed at RHIC. 



VI. CONCLUSIONS 

MACE shows that if sound waves are created by a 
jet propagating through Quark-Gluon-Plasma, two- and 
three-particle correlation functions will show the ex- 
pected double- and fourfold-peaked structures. The cone 
is not washed out by collective flow. 

The peaks will, though, not be at the expected "perfect 
mach positions" A(fi,2 = 7r±cos _1 cs, but shifted to 
positions expected for lower sound velocities. This shift 
of the máxima is strong enough (5A.ip max ss 0.2) to feign 
the presence of Hadron gas in favour of the (actually 
present) Quark-Gluon-Plasma. 

The presence of a double-peaked structure (for two- 
particle correlations) means that no wide-peak structure 
arises in MACE. That means that it is not created by an 
enormous shift of the conical angles by flow. Non-linear 
effeets in the cone which might oceur (spatially) near the 
jet may be a better explanation. If the created shape is 
not a true cone, but a cone with a round top, partióle 
emission will be enhanced in backward región with re- 
spect to the results from MACE. Also, a slowing down of 
the jet will change the correlation signal to enhance the 
backward región. 

Event-by-event analysis might give better insight into 
the real speed of sound. Correlation functions from mid- 
rapidity jets that are symmetric around n show their 
peaks at the expected positions. But single-event trig- 
gers may be useless due to event-by-event fluctuations. 
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